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Abstract. We prove existence and uniqueness of stochastic equihbria in a class of 
incomplete continuous-time financial environments where the market participants 
are exponential utility maximizers with heterogeneous risk-aversion coefficients and 
general Markovian random endowments. The incompleteness featured in our set- 
ting - the source of which can be thought of as a credit event or a catastrophe - 
is genuine in the sense that not only the prices, but also the family of replicable 
claims itself is determined as a part of the equilibrium. Consequently, equilibrium 
allocations are not necessarily Pareto optimal and the related representative-agent 
techniques cannot be used. Instead, we follow a novel route based on new sta- 
bility results for a class of semilinear partial differential equations related to the 
Hamilton-Jacobi-Bellman equation for the agents' utility-maximization problems. 
This approach leads to a reformulation of the problem where the Banach fixed 
point theorem can be used not only to show existence and uniqueness, but also to 
provide a simple and efficient numerical procedure for its computation. 

1. Introduction 

Market incompleteness and equilibria. The central theme of this paper is a study of an 
equilibrium problem in an incomplete continuous-time stochastic setting. In contrast with the 
complete case where significant advances have been made in continuous time (see, e.g., [11171 [^[TUl 
[TTJ [ini [T71 [THl [31] as well as Chapter 4. of \i9l) the incomplete-market literature is lagging behind. 
Even among the few incomplete markets studied so far (see [2] or [14] , for example) ideas related to 
market completeness, typically through the representative-agent approach, are used. To the best 
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of our knowledge, the present paper is the only one in continuous time where a fully-incomplete 
market structure, in the sense that both the prices and the set of replicable claims (the marketed 
subspace) are determined as a part of the equilibrium, is analyzed and existence of equilibria is 
established. 

The difficulty with incomplete markets is that Pareto optimality, which is commonly exploited to 
establish existence of equilibrium in a complete market, is not guaranteed anymore. Our approach 
rests upon the notion of stability of demand - a continuity property of the optimal response (optimal 
portfolio) when viewed as a function of market dynamics (the market-price-of-risk process) . The 
leitmotif behind such an analysis is the following: if the aggregate demand can be shown to possess 
good continuity properties in an appropriate topological setting, a fixed-point-type theorem can be 
used to guarantee existence and (with some luck) uniqueness of an equilibrium market dynamics. 

Stability in Holder spaces. The technical bulk of the present manuscript is devoted to the 
stability of the optimal investment strategy under small functional perturbations of the market- 
price-of-risk coefficient. Problems similar to ours have attracted some attention recently - see, for 
example, [H O [HI [201 [23l [24] . There are at least three - at first glance very different - reasons why 
such problems are important. 

First, from the statistical and financial point of view, it is important to understand how misspec- 
ification or misestimation of the market dynamics coefficients affects the optimal trading strategies 
of agents who take the coefficient estimates at face value. Equivalently, one can wonder what 
kinds of statistical procedures for the estimation of those coefficients yield the most stability in 
implementation. For a deeper discussion of this point, we refer the reader to [24|. 

Second, in agreement with the classical methodology of the theory of partial differential equa- 
tions, and applied mathematics in general, the following three aspects of every new problem are 
typically studied: existence, uniqueness and sensitivity of the solution with respect to changes of 
the problem's input parameters. These criteria are generally known as Hadamard's well-posedness 
requirements (see [13j). We view model specification as one of the most important input data in the 
utility-maximization problem, and understand the stability with respect to it as one of Hadamard's 
requirements. 

The last reason - by far the most important one for the present paper - is to shed more light 
on the intimate relationship between the notions of stability and that of competitive equilibrium 
in incomplete markets. Indeed, rationality in the presence of individual preferences can be viewed 
as the "law of motion" of financial agents. In aggregate over several investors, the optimal trading 
strategies can be interpreted as the agregate demand for the assets in question. This aggregate de- 
mand, via the principles of market clearing, ultimately determines the shape of the price-dynamics 
of the financial assets. 

Unlike other stability-related results mentioned above, the present paper approaches the problem 
from a different point of view. While all the previous strategies involved probabilistic and convex- 
analytic techniques, we tackle the problem from a pure PDE perspective. In this way we are 
able to draw much stronger and more precise conclusions about the behavior of the function 
which maps the market-price-of-risk process into the optimal portfolio. Our main stability result 
is that this map is locally Lipschitz continuous when the inputs and the output are placed in 
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anisotropic Holder spaces on [0, T] x M. and we give explicit estimates on the Lipschitz constant 
(see Theorem 12. 6|) . The price we pay for the PDE approach is not as dire as one may think when 
placed in the equilibrium setting - the necessary Markovian assumption only restricts the form of 
the agents' random endowment to functions of state variables. The Markovian structure of the 
resulting equilibrium asset dynamics is, on the other hand, not a restriction. It can be viewed as a 
strenghtened version of an equilibrium existence result: not only do we show that an equilibrium 
exists, we also show that a Markovian equilibrium exists. 

In order to illustrate our techniques with the minimal amount of distraction, we choose what 
can be termed as the simplest nontrivial incomplete market model. Indeed, our market consists 
of a single risky asset (and a unit riskless asset) whose dynamics depend on a single Brownian 
motion and a one-jump Poisson process. All the incompleteness in the market comes from this 
uninsurable and unpredictable jump. Every agent is an exponential utility maximizer who also 
receives a random endowment at the end of the time horizon and this random endowment is allowed 
to depend on all sources of uncertainty in the market, including the indicator of the unpredictable 
jump. As a form of a normalization, we assume that the volatility of the risky asset is constant, 
but the drift can depend on time and the current value of both state processes - the Brownian 
motion and the indicator of the unpredictable jump. It is clear that our setting can be generalized 
in many different directions and that our estimates are far from optimal. Our goal was to provide 
a proof of concept for the powerful PDE techniques outlined in the body of the paper. We strive 
to keep the presentation as simple as possible and as accessible as possible to a reader who is 
not a specialist in PDE. For that reason, we use only elementary techniques and results in the 
Schauder-type theory of semilinear parabolic PDE as outlined, for example, in [21^ . In the same 
spirit, we do not pursue any connections with BSDE. 

Existence and uniquencess of the equilibrium price dynamics. Once the local Lipschitz 
continuity of the agents' demand functions is established, we turn to the market-clearing conditions 
and show that they can be rephrased in terms of a simple fixed-point problem for a continuous 
map on an anisotropic Holder space. Surprisingly, under a "smallness" condition, a restriction of 
this map becomes a contraction and therefore, both existence and uniqueness of equilibrium can 
be guaranteed. Furthermore, an efficient and easy-to-implement numerical technique - based on 
iteration - emerges naturally. Results of this type are very rare in equilibrium theory; equilibria 
are typically not unique and their multiplicity is often one of the major objections to the use 
of equilibrium modelling in practice. Also, the computational methods (even in the fully-finite- 
dimensional case) are often quite involved and based on inefficient and relatively hard-to-implement 
procedures based on Sperner's lemma and the related Scarf's algorithm (see [S^ and |29]V 

The structure of the paper. Section 2. desribes the model and provides statements and some 
proofs of our main results. Section 3. is devoted to the proof of the central stability results, while 
Appendix A contains the pertinent information about Holder spaces. 

2. The financial environment and market clearing 

We define our financial environment by specifying the structure of the three main ingredients: 
the information structure, agents' preferences and the completeness constraints. 
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2.1. The information structure. Let T > be a real number which we interpret as the time 
horizon, and let {Bt}te[o,T] be a standard Brownian motion defined on a complete probability space 
{fl, T, P). We assume that, additionally, the same probability space accommodates an independent, 
exponentially distributed random variable r with parameter /x > 0, where the corresponding 
counting process {]^t\t^^,T\ is given by 

iVt = l{^<t}, [0,T]. 

Let {J"t}4g[o be the right-continuous augmentation of the filtration \Tf V .^t^}te[o,T]) where 
Tf = a({BJ,<t) and = a({iVJ,<t) for t € [0,T]. 

A more detailed interpretation of the roles of the two filtrations will be given later. For now, 
let us just mention that the tr-algebra J* models the total publicly-available information at time 
t and that r models an event with a significant impact on the market. Our model does not deal 
with asymmetric-information situations - all the agents have access to the same information, and 
no information is hidden from any of them. As we shall explain below, it is the way in which the 
information trickles into the dynamics of traded assets that makes the market incomplete. 

2.2. Completeness constraints. The final goal of our analysis is the determination of the form 
of an equilibrium asset-price process. Without exogenously-imposed constraints, there is nothing 
that will prevent the agents from "opening" as many markets as possible, and, eventually, building 
a market structure that will be able to replicate any uncertain pay-off. In other words, without 
constraints - and bar pathologies - all equilibrium markets are necessarily complete. The situation 
we are modelling, however, calls for a degree of incompleteness: we envision the situation in which 
the random variable r marks an event with the property that the time-scale at which the market 
adjusts to T is larger than the time-scale on which it reacts to the other information. Consequently, 
as we shall see, its nature will be such that no perfect insurance against its effects can be bought 
or sold. 

More generally, market environments where a portion of the information fiow affects prices 
faster than the rest are important examples of natural constraints that preclude completeness; 
we term them fast-and-slow-information market environments. In fact, many widely used 

incomplete market models (virtually all models where the incompleteness stems from more "sources 
of uncertainty" than available risky assets) can be viewed as equilibria in fast-and-slow-information 
market environments. The situation described in the present paper corresponds to, arguably, the 
simplest continuous-time fast-and-slow-information market environment. More complicated (and 
realistic) environments can be constructed and analyzed, but we opt for one which allows us to 
showcase our methods without unnecessary confusion. 

The defining properties of r lead to the following class of possible asset-price dynamics: let A be 
the family of all { J"t}tg[o^T]-predictable processes {)^t}te[o,T] such that \Xu\ du < oo, a.s., and 
let <S = |{'5't^'}te[o,T] : A € a| be the family of all Ito-processes {-St^'jtgio.T] with the dynamics 
given by 

dSl^^ =Xtdt + dBt, S^Q^^ = 0, (2.1) 
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as A ranges over A. Additionally, let A'^'^ be the subset of A consiting of all uniformly bounded A. 
The corresponding set of price processes is denoted by S^'^. We always (implicitly) assume that a 
risk-free numeraire asset S = 1 accompanies each S'-^\ 

Remark 2.1. 

(1) It is important to note that there is no loss of generality in assuming that the volatility 
coefficient is identically equal to 1 or that "arithmetic" dynamics for 5^'^^ are chosen instead 
the more common "geometric" one. Indeed, the only property of the asset-price dynamics 
important for the determination of the equilibrium is the set of replicable random variables 
it produces, and this depends on drift and volatility only through their quotient A (the 
so-called market price of risk). In fact, our formulation of the equilibrium simply cannot 
distinguish between the two. This would be the case in the more complicated (and more 
realistic) model where the agents' random endowments depend on the asset-price itself, 
i.e., where the agents hold financial derivatives of the asset. Similarly, the postulated 
availability of a trivial numeraire asset 5* involves a minimal loss of generality. Indeed, 
the agents' utilities depend solely on the terminal wealth and do not exhibit any time- 
impatience characteristics. In fact, without introducing consumption into the model, it is 
impossible to disentangle the two solely on the basis of equilibrium analysis. 

(2) The two terms on the right-hand side of (12. ip contribute to the fluctuations in 5^'^^ on 
different scales; the order of magnitude of the term Xdt is dt while the order of (the 
absolute value of) the term dBt is Vdi. Equivalently, the first term takes on average 
l/\/dt times longer to produce the same (local) efect as the second one. In this sense, the 
fact that the new information regarding r comes in only through the first, slower, term 
corresponds exactly to the requirement that the occurrence of r be absorbed into the asset 
dynamics on a slower scale than the shocks produced by B. 

(3) One can view (|2.ip as a restriction on the set of stochastic processes that may serve as 
allowable market dynamics. For a different specification of the characteristics of the mar- 
ket environment (constraints on the information flow, number of assets, etc), we would 
get different families of processes. Therefore, one could give the following, abstract, defi- 
nition: a completeness constraint is simply a family S of (possibly multi-dimensional) 
semimartingales . 

(4) While no dividend structure has been mentioned so far, there is a way of incorporating 
dividends in the present model, and, perhaps, the best way to describe it is to compare it 
to a particular incomplete setting - the one with short-lived securities - in the language of 
[5S] (see also how it leads to existence results in discrete time as in [5S], Proposition 25.1, p. 
255). One can think of the asset S*'^^ as an aggregation of short-lived assets with inception 
t and maturity t + dt, each of which pays a dividend Sj:^^^ and costs Sf^\ Equivalently, 
the agent enters into a bet a unit of which pays dBf and costs —Xfdt. In this way, we can 
think of dBt as the local dividend and —At as its price. It should be noted that it is easy 
to incorporate any (reasonably behaved) martingale Mt of the form Mt — Jq CTu dB^ as the 
divided process - one has to replace S^^^ by Uu dS^\ Even if M is not a martingale, we 
can simply add its drift to — A and turn it into one. This procedure also allows us to give 
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a loose interpretation of the formation of volatility in the market. While the market price 
of risk A is determined in the equilibrium from the agent's primitives and the structural 
properties of the dividends, the volatility is determined by their quantitative properties. 
Consequently, the unit-volatility assumption we impose on S^^^ can be reinterpreted as a 
simple normalization of the dividends. 

2.3. The agents and their preferences. We assume there is a finite number / G N of agents, 
all of whom actively participate in trading in all available assets. The preference structure of each 
one of them is determined by the following: 

(I) the utility functions: each agent is an exponential-utility maximizer with 

W{x) = -exp(-7*a;),a; € R, 7' > 0. 

(II) the random endowment: f = g'^{BT, Nt), with bounded : M x {0, 1} M, where further 
regularity conditions are to be specified. 
As usual in the Alt-von Neumann-Morgenstern expected utility paradigm, agent i prefers the 
random variable Xi E Ft to a random variable X2 G Ft if and only if E[J7* [Xi +f )] > E[J7* {X2 + 
f )], where we set E[F] = -co whenever E[y+] = E[F^] = +00. 

Remark 2.2. 

(1) If one thinks about B and N as factors, (II) above states that the agents' random endow- 
ment depends on all factors driving the public information. In particular, the (conditional) 
distribution of the random endowment may change abruptly and considerably at time r. 
One of the possible financial interpretations of the situation is that all agents hold (long 
or short) positions in assets whose pay-offs are affected by the occurrence of r (default- 
sensitive derivatives, callable bonds, disaster-sensitive investments, etc.). 

(2) The exponential nature of the utilities allows us to partially remove the assumption that all 
agents use the same (subjective) probability to compute the expected utility of a particular 
position. Indeed, using the identity 

E'" [- eM-l\X + r))] = E[- exp(-7''(X + T))], 

where = 5* — i log(^), we can easily "absorb" different subjective probabilities into 
the form of the random endowment. Care must be taken, though, to ensure that the 
appropriate integrability conditions are met. Also, it should be noted that such a change 
of measure can lead to loss of the Markovian structure of the ingredients. 

Let us now focus on the pertinent case when the set of tradable assets consists of a single risky 
asset given by S^^'' e S^'^, for some bounded market-price-of-risk process A € A'''^ and the trivial 
numeraire asset S (as described just below equation (|2.ip . above). Agent i chooses a dynamic 
self-financing portfolio process in the appropriate admissibility class (to be specified shortly) so as 
to maximize the expected utility of terminal wealth: 

WUq dS\^^ -f £■')] max . (2.2) 

Here, the value ttj of the one-dimensional process {T^tltel^.T] denotes the number of shares of the 
risky asset in the portfolio. We do not explicitly mention the number of shares pt of the riskless 
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asset, since, thanks to the self-financing condition, it necessarily satisfies 

Pt^ f\udS\^^ -^tS[^\ (2.3) 

Our class of admissible portfolio processes aims to be just restrictive enough to rule out doubling 
strategies, yet large enough to contain the maximizer of the utility-maximization problem p.2jl . 
We would like to emphasize that, due to the regularity of some of the ingredients, one does not need 
the sophistication typically encountered in general semimartingale models (see, e.g., the classes 0^, 
I = 1, 2,3,4 in [8J or the notion of permissibility in [26]). Instead, we sacrifice a small amount of 
generality for a large gain in simplicity by proceeding as follows: the admissibility class A consists 
of all {J"t}jg[o,T]-predictable processes {T^t}t£[o,T] such that E[exp(jJ'(i -|-e)7r^ du)] < oo, for some 
£ > 0. Note that, by Holder's inequality, the integral jj" 7r„ dS^^ is well-defined for each -k £ A 
and A G A. 

Remark 2.3. As we shall see later, the optimization problem (j2.2|) admits a, dt x dP-unique max- 
imizer TT*^^''* of (|2.2p in the class A for each A G A*^^ and each agent i — 1, . . . , /. The portfolio 
process {''^i^^'^}t£[o.T] is interpreted as an optimal response of agent i to the market dynamics 
induced by A; the mapping {^t}te[o,T] ^ {^t^'''*}tG[o.T] pl^^ys a central role in our analysis. 

2.4. Market-clearing conditions and the main result. Having introduced all the ingredients, 
we turn to our main problem. A central economic principle is that the prevailing market dynamics 
must have the following fundamental property: the demand and supply for each tradable asset 
must match, i.e., all markets must clear. A precise definition follows: 

Definition 2.4. A process S''^'^^ e S (or, equivalently, its market-price-of-risk process A e A) is 
said to be an equilibrium price dynamics (market price of risk) if there exist processes 
{7rj'^''''}fg[o.T] €z A, 1 = 1,...,/ such that 

(1) (rationality) for all vr G ^ and all for alH = 1, . . . , /, 

nU'iJo T^u^'^ dSi^^ + £■')] > E[W{J^ 7r„ dsi'''' + £% and 

(2) (market clearing) J^Li 4^'''* = 0, for all t e [0,T], a.s. 

Remark 2.5. The reader will immediately realize that we only require market clearing for the risky 
asset S'''^\ The self-financing condition (|2.3p implies, however, that the market for the riskless 
asset S will clear, in that case, as well. 

The main result of the present paper - which asserts the existence, uniqueness and efficient 
computability of equilibrium price dynamics - is summarized in Theorem 12.61 below (we direct the 
reader to Appendix |A] for details on Holder spaces and the related notation) : 

Theorem 2.6 (Existence, Uniqueness and Computability of Equilibria) . Consider the setup given 
in subsections \2.1{ \ 2.2\ and \2.3[ and assume, additionally, that g^{-, n) € C^^"(R), for i = 1, . . . , I , 
n = 0,1 and some a £ (0,1]. Then there exists a constant Tq > (which may depend on 
a^M, {7'}i=i,...J o-nd {W\2+a)i=i,...,i) such that for T < Tq, 

(1) the equilibrium price dynamics exist and the corresponding market-price-of-risk process is 

unique in the class of all processes {^t}t&[o,T] such that At — X{t, Bt, Nt-), t G [0,T], for 

some function A G C'^+°'{[Q,T] x M x {0, 1}). 
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(2) the function X is efficiently computable, i.e., there exists a subset D C C^^"([0,r] x M x 
{0, 1}) with Cz D and a contraction II . D ^ D such that X is the fixed point of II. 



(1) The most restrictive feature of Theorem 12.61 is the "smallness" condition we place on the 
size T of the time horizon. It is not evident that it can be dealt with by "pasting" equilibria 
together as the upper bound on T depends on the primitives of the model. In particular, 
extending an equilibrium from T to T' > T would entail using the value of the original 
equilibrium at time T as the terminal condition of the equilibrium on [T, T'] . An iteration 
of this procedure may very well lead to a bounded sequence. 

We do conjecture, however, that this assumption can be relaxed if one only wants to 
establish the existence of the equilibrium, but cannot be relaxed if the additional benefits of 
uniqueness and efficient computability are also desired. The resolution of those important 
questions is the content of our future work and seem to require quite different and more 
sophisticated techniques. Moreover, preliminary numerical experiments suggest that the 
for a large class of realistic parameter values, the smallness constraint does not seem to be 
binding and the iteration procedure converges. 

(2) It is, perhaps, instructive to interpret the above uniqueness result in the light of part (4) 
of Remark 12.11 What is truly unique is the equilibrium price of the local gamble of size 
dWt . It follows that each replicable derivative security C admits a unique equilibrium price 
process S. The difference with the complete models is that the class of replicable contingent 
claims is determined as a part of the equilibrium and admits no simple description a priori. 

The proof of Theorem 12.61 is based on the following stability estimate which, as discussed in 
the introduction, is of interest in its own right. We devote the entire section [3] to its proof. 
Since it deals with a single-agent's optimization problem, we omit the agent index i from its 
statement. Also, in the interest of readability, we introduce some additional notation: for a 
function u : [0, T] x R x {0, 1} — > M, we define the n-difference it„ as 



For such functions, restrictions •, 0) and •, 1) are called the n — {)- and the n — 1-slices of 
u. We say that a measurable function / : [0, T] x M x {0, 1} — > M is a Markov representative of 
the process {l^t}te[o,T] if 



It is convenient to identify a process admitting a Markov representative with the Markov repre- 
sentative itself. In fact, we will do so from this point on with little or no explicit mention. 

Theorem 2.8 (Stability of the Optimal Portfolio). For X G C"([0,r] x M x {0,1}) and g € 
C2+"(M X {0, 1}) there exists a unique function u^^^ G C^+°'{[0,T] x M x {0, 1}) which solves the 
following Cauchy problem for a semilinear PDDE (partial diferential- difference equation): 



Remark 2.7. 




Yt^ f{t,Bt,Nt-), for alHG [0,T], a.s. 




on [0,T) X M X {0, 1}, 



(2.4) 
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The stochastic process {'^['^''}te[o,T] with the Markov representative n^^^^ = ^-^ ~ Ux , i-e., the 
process given by 

ni^^ ^^X{t,Bt,Nt.)-ui^\t,Bt,Nt^), t e [0,T], (2.5) 

is in A and attains the maximum in the utility-maximization problem (|2.2p . 

The maps A i— > tt*^'**-', A i— > Ux^\ from C"([0,T] x R x {0,1}) to itself, are locally Lipschitz 
continuous. More precisely, there exists a constant C = C{'y,fi,a) > such that for each R> 
and Ai,A2 G C"{[0,T] x M x {0,1}) with \X2\a < R we have 

\ui^'^ - ui^'^U < L{R)\X, ^ and|7r(^i)-7r(^^)|„< (i + i(i?)) IA1-A2U, (2.6) 

where 

L{R) = CT(i+")/(2+")e^^"^^""°"""'"^^" + (1 + T)(l + R^))'^'". (2.7) 

With Theorem 12. 81 at our disposal, the proof of Theorem [53] becomes straightforward. We start 
with a characterization which follows directly from Definition 12.41 and equation (|2.5p of Theorem 
12.81 When the index i is added to u'^^ (or its derivatives), as in u^^'^'^ (or ui^^'^, etc.), we are 
referring to the solution to (j2.4p with the agent-dependent terminal condition g — and the 
risk-aversion parameter 7 = 7*. 

Lemma 2.9. A process {Xt\te[vi.T] with Markov representative X G C"([0,T] x M x {0,1}) is an 
equilibrium market-price- of -risk if and only if it is a fixed point of the operator 11 : C"([0, T] x M x 
{0, 1}) ^ C"([0, T] X M X {0, 1}) defined by 

I 

n(A) = 7E"i'^'' 

i=l 

where ^ ^ (Y.\=i^) ■ 

In the sequel, let i?o = ^ 1^2;°'' (where ui^"*'* corresponds to A = 0) be computed for 

T = 1. Note that i?o dominates the value of the same expression we get when we set T < 1. 

Lemma 2.10. There exists a constant Ti — Ti (a, /i, {7*}^=!, {|f;'|2+a}i=i, ...,/) > 0, such that 
ifT<Ti, n maps the ball 

B,(i?o)={AeC"([0,r]xMx{0,l}) : |AU < i?o}, 

of radius Rq in C"([0, T] x R x {0, 1}) into itself 

Proof. The first inequality in (|2.6p of Theorem 12.81 implies that there exists a non-decreasing 
function F : [0, 00) — [0, 00) such that if T < 1 we have 

l4^^1a < ki^^la + ri/(2+")^^(|AU), and so |n(A)U < ii?o + i/Ti/(2+")i^(|A|„), 

for A e C"([0,r] X R X {0,1}) and alii = 1,...,/. Note that F, being derived from L(-) of 
(|2.7p . depends on a, ^, 7', |5*|2+q, i = l,---,-^ but not on T or A for T < 1. Therefore, for 
T < Ti = 1 A ( amfl ) )^"^°' have |n(A)|„ < Rq, whenever |A|a < Ro- □ 
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Lemma 2.11. There exists a constant Tq = To^a, fi, {\g^\2+a}i=i....j) > such that 

if T < To, the mapping 11 is a contraction from Ba{Ro) into itself. 

Proof. By Lemma [2. 101 the mapping 11 maps Ba{Ra) into itself, as long as T < Ti. Therefore, we 
can use the Lipschitz estimate (|2.6p of Theorem 12.81 and reasoning similar to the one in the proof 
of Lemma 12.101 to conclude that there exists Tq < Ti such that for T < Tq, we have 



Proof of Theorem ] 2. 6\ The path to the proof is paved by Lemmas 12.91 12. 101 and 12.111 Indeed, we 
simply apply the Banach fixed point theorem to the mapping 11 on the complete metric space 
Ba{Ro) and note that, by the first part of Theorem[5^ n(A) € Ci+"([0, T] x R x {0, 1}), as soon 



The purpose of the present section is to prove Theorem 12. 81 We break its statement into several 
smaller, more manageable results, and proceed to discuss them one by one. 

3.1. The HJB equation: existence and verification. The first assertion of Theorem 12.81 
is that a the PDDE in (j2.4[) admits a regular solution and that it can be used to construct a 
Markov representative for the optimal portfolio in the utility maximization problem (j2.2p . Similar 
characterizations are ubiquitous in stochastic control in the case of exponential utility, and are too 
numerous to list; we simply instruct the reader to consult [6 and references therein. We note that 
a PDDE is a special case of a partial integro-differential equation, but we choose not to use that 
term because of the simplicity of the integral component. Equivalently, it could have been stated 
as a system of one linear and one semilinear PDE. 

Appendix |A] should be consulted for notation not explicitly introduced in the main body of the 
paper. However, for the reader's convenience, we state the following convention both here and in 
subsection IA.2I of Appendix 1X1 

Convention 1. 

(1) The variables fx, 7 and a G (0, 1) are considered "global" and will not change throughout the 
paper. Any function of the global variables (and global variables only) if called a universal 
constant. 

(2) The notation a < b means that there exists a universal constant C > such that a < Cb. 
Such a constant may change from line to line. 



Our analysis starts with a simple linear existence result in the spirit of Schauder's theory: 
Lemma 3.1. For h,ae C"{[0,T] x R) and g e C2+"(R) the Cauchy problem 



|n(Ai)-n(A2)U<i|Ai-A2| 



for Ai,A2 e Ba{Ro)- 



□ 



as A € C"([0,r] X R X {0,1}). 



□ 



3. Stability of the optimal portfolio: a proof of Theorem 12.81 




(3.1) 



(3.2) 
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Proof. The existence of a unique solution to dSTj) in C^+"([0,r] x R) is well-known (see, for 
example (21], Theorem 9.2.3, p. 140). To get p.2p . we use the fact that if u is the unique solution 
to ([37T|) . then u{t,x) = e"'^(^~*)u(i, x) solves 

{0 = Ut + hiixx — 0u + hu^ + a on fO, T) x R 
(3.3) 
^(r,-) = .9(-), 

where a{t,x) = e"''('^"*)a(t, a;). It remains to note that the constant function w(t,x) = \g\Q + 
^^l^'lo = \g\o + is a subsolution, and that its negative — w is a supersolution of (|3.3|) : the 

maximum principle (see |21| . p. 105, Theorem 8.1.2) implies 

-w{t, x) < u{t, w) < w{t, x), {t, x) e [0, T] X R. (3.4) 

□ 

When a non-linear term is added, a similar result can be obtained with a bit more work. The 
argument is based in part on the following lemma: 

Lemma 3.2. Let {a;„}„gNo be a sequence of nonnegative real numbers with the property that 

xn+i <A + Bx^, (3.5) 

for some constants A, B > and < a < 1. Then there exists a constant g{a) > 1, independent 
of A, B and {xn}neN> such that 

limsupxn < 5(a) max(^, B^/*^^^"^). (3 5) 

n 

In particular, the sequence {xn}neN is bounded. 

Proof. If B — 0, the inequality in (|3.6p clearly holds. When A = and B > 0, we have Xn < 
^i+a+a^+ . +a'-i^;""^ which implies (ESI directly. 

Focusing on the case A,B > 0, we set B — B/A^^", and let G > 1 be the unique positive 
solution to G = 1 + BG". The scaled sequence ?/„ = x„/(AG), n e N satisfies 

Vn+i = < 1/G + BG"-\ynr = K + (1 - K)y^, (3.7) 

where k = G"^ = 1 - BG"^^ . Therefore, if ?/„ < 1, then yn+k < 1, for all fc e N. On the other 
hand, suppose that y„ > 1, for all n € Nq. Then, for n € No, we have yn+i < Vn^ so, y„ < j/g". 
Consequently, lim sup„_j.Q^ y„ < 1, i.e., limsup„a;„ < AG. 

It remains to show that AG is bounded from above by the expression on the right-hand side of 
(|3.6p . Let g{a) be the unique positive solution of g{a) = 1 + .9(0:)", so that g{a) > 1. If i? > 1, 
then 

1 + B(,g(a)Bi/(i-"))" = 1 + {g{a) - 1)51/(1-") < g(a)B^^^^~°'\ 
Therefore, the monotonicity of the function x x — I — Bx°' implies that G < g{a)B^^'^^^°'\ 
When _B < 1, we have 1 + Bg{a)°' < g{a), so G < g{a). Therefore, 

AG < Ag{a) max(l, ^ g(^a) max(A, B^/^i-"'). 

□ 
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The existence of a solution to semilinear equations of the type treated in Proposition 13.31 is 
well-known (see, for example, [3]). However, the C^+"-regularity and the |-|q and [•]2+q estimates 
need additional work in this proof. The reader will note that we still go through the steps of 
the existence argument; the reason is that the Holder regularity is established through the | • |o- 
convergent sequence which is originally constructed to show the existence of a solution. 

Proposition 3.3. Consider a semilinear Cauchy problem of the form 

( = ut + hu^,x + hu^ + a - 5e''" on [0, T) x M, 
y u[T, ■) g, 

where g G C2+"(M), h £ C"([0,T] x R), a,b e C"([0,r] x M), 7 > and a{t,x) > 0, b{t,x) > 0, 
for all {t,x) e [0,T] x R. Then i3.8\) admits a unique solution in C^"'""([0,T] x M). Furthermore, 
the following bounds hold for all {t,x) G [0,T] x R; 



_i log(e''lf lo + jT\b\o) < u{t, x) < \g\o + T\a\o, (3.9) 
[uh+o. ^ [gh+o. + Ha + [6]„e'^(lsl"+^l''l«) + (1 + |5|i+"e^(i+i")(^l-l»+l9l«) + \h\l+")\u\o. (3.10) 



Proof. We start by setting F{t,x,y) = a{t,x) — a;) exp(7?/), and note that, by Lemma I A. 6 1 
F{-, ■,w{-, •)) is in C"([0, T] xR), whenever w is. Therefore, by LemmaO for each w e C"([0, T] x 
R), there exists a unique C^+"([0,r] x R)-solution u of the Cauchy problem 

(0 = ut + ^Uxx + hux + Fi-,-,w{-,-)) on [0,T) x R, 

\u{T,-)^g{-). 

Consequently, due to LemmaE] the operator G : C"([0, T] x R) ^ C'^+°'{[Q, T] x R) which assigns 
to w e C"([0,T] X R) the unique solution of ([XTT]) is well-defined. For ui,U2 e C"([0,r] x R), 
the function k — Qui — Gu2 satisfies 

J = Kt + ^K^.^ + + 6 on [0, T) x R, 

|«;(T,.) = 0. ^^'^^^ 
where S{t,x) — F{t,x,ui{t,x)) — F{t,x,U2{t,x)). The inequality p.2p of Lemma 13. II implies that 

|k|(/3) = \GU2 - e^Ul|(/3) < }\S\^13). (3.13) 

On the other hand, using inequality (jA.10|) of Lemma [A. 61 we obtain 

\Syj < /3e-^(^-*) |6(t,x)|7e^r"ilvM |^^(^^^) _ u,{t,x)\,, (3.14) 

where we remind the reader that \u\ = sup^^ x)e[o t]xr^(^i^)- 

If u belongs to the range TZQ of G, then the maximum principle applied to equation (|3.1ip and 
the non-negativity of b imply that 

M < r|a|o + Islo- (3.15) 
Inequalities p.l3|) . (|3.14p and (|3.15|) show that for ui,U2 € TZG 

\gu2 - Gu^y^ < i|6|o7eT(^l''l«+l^'l°)|u2 - (3.16) 
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It follows that the mapping Q is a, \ ■ |(^)-contraction for (3 > large enough. For such (3, the 
sequence where ui € C"([0,T] x R) and Un+i = Gun, n Cz N, converges towards some 

u G C([0,r] X M) in I • l(^), and. therefore, also in | • |o. Our next task is to show that u is not only 
in C([0, T] X R), but also in C"([0, T] x R) and that it is, indeed, a fixed point of g. 

Let u e TZg and set u = Qu. Then, the inequality (|3.15[) is satsified by u as well. Moreover, 
applying the maximum principle once again to equation (|3.11|) , we get that 

-Islo - r|6|oeT(^l"l°+l9l°) < u{t,x). (3.17) 

Combining p.l7p with inequality p.lSp for u, we obtain 

|w|o < |.9|o + Tmax{|a|o, |6|oe^(^l''l''+l9l")}. (3.18) 

which provides a uniform | • |o-bound on all elements of the sequence {un}„gN- In the remainder of 
the proof, D will denote a generic constant which may depend on a,T,j, \h\a, \a\a, \b\a or |(7|2+a, 
but is independent of n and ui and may change from occurrence to occurrence. Thanks to the 
uniform bound on |m„|o established in ()3.18|) and recalling that Un+i = Gun G C^+"([0,r] x R), 
the first inequality in (jA.Sjl of Theorem IA.3I yields 

[un+iU < D[u^+,]^i^^+"\ (3.19) 



On the other hand, CoroUarv I A. 5 1 implies that 

K+i]2+a <^([5]2+a + (l + |/j|'+")|«„+i|o + [a-^'e''""]„), (3.20) 

while by inequality (IATO)) of Lemma [AJl we have [e'^'"]a < 7e'^r«'l for any w £ C"([0, T] x R). 
So, using once more the uniform bound (|3.18|) on |u„|o, we get 

[Un+lh+a < D{1 + |5|a|eT""|„) < D{1 + K]„). 

Therefore, by (|3.19p . we have 

< D{1 + K]a)"/(2+"), for n e N. 

This fact implies that the sequence Xn = [un]!^^"''^" satisfies p.5[l with A — B — d(2+q)/q ^-j-^g^^ 
by Lemma [3.21 

SUp[M„]a < OO. 
n 

It is not difficult to see that the | • |o-closure of a | • ja-bounded subset of C([0, T] x R) is, in fact, 
a subset of C"([0,T] x R). Consequently, the | • jg-Umit u of {u„}„eN is in C"([0,r] x R), and in 
particular, Qii is well-defined. Inequality p.l6p yields 

IGu - M„+i|o = \Gu - GUnlo < D\u - Un\o, 

SO that Gu — lim„ u„ — ii, i.e., ii solves ()3.8p . Moreover, any solution u must satisfy (I3.15P and 
the relation ()3.16|) guarantees that u is unique in the class C"([0,T] x R). 

To establish the bounds in p.9j) we pick two continuous functions a, b : [0, T] — >■ R and consider 
the function 

w{t, x) = A{t) log \jB{t) + exp(-7G') 

7 L 
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where A(t) — a{u) du, B{t) — b{u)e'^'^^^^ du > 0, and G is an arbitrary constant. As the 
reader can easily check, we have 

=Wt + ^Wxx + h(t, x)wx + a(t, a;) — h{t, x)e''^ 

- {a{t,x) - d{t)) - (b{t) - 6(i,a;))eT"'(*'") 
w{T,-)^G. 

For different choices of functions a, b and the constant G, w will be either a sub- or a supersolution 
of (|3.8p . Indeed, for G — \g\Q, d{t) = |a|o and b{t) = 0, w becomes a supersolution, yielding the 
upper bound in (|3.9|) . Similarly, for G = — Ifflo, d{t) — and 6(t) = |6|o, ui is a subsolution, and 
so, the lower bound in (|3.9p holds, too. 

The last item on our list is the [•]2+Q-bound p.lOp . Estimate (jA.7|) . along with Lemma [A. 61 and 
the just established (|3.9p , yields 

M2+C ^ [gh+c. + [a]« + [te^la + (1 + |/i|'+")|u|o 

^ [gh+c + Ha + [felae-'r^T + |6|o[eT"]„ + (1 + \h\l+'')\u\o 

^ [gh+a + [a]a + [5]„e'^(^l''l«+lffl°) + |6|o7e^^^'"'°+'^'°^Ma + (1 + |/i|'+")|m|o. 

The inequality (|3.10p now follows directly from Corollary lA. 21 □ 

Finally, we integrate Lemma [XT] and Proposition 13.31 and discuss the difficulties in the standard 
verification argument. 

Proposition 3.4. For X e C"([0, T] x M x {0, 1}) and g e C'^+"{R x {0, 1}) there exists a unique 
function u(^) £ C'^+"{[0,T] X M X {0, 1}) which solves the Cauchy problem Moreover, the 

stochastic process {T^i^'^ytelo.T] with Markov representative n^^^ = — ui^"^ is in A and attains 
the maximum in the utility maximization problem (|2.2p . 

Proof. We note, first, that the n = 1-slice of the equation ()2.4|) has the form (|3.1j) with h{t, x) = 
-X{t,x) e C"([0,T] X R), a{t,x) = ^X'^(t,x) G C"([0,T] x R). By Lemma [SH it admits a 
unique C2+"([0,T] x Resolution Therefore, the n = 0-slice is of the form with 

h{t,x) = —X{t,x), a{t,x) — i(iA^(t, a;) + /i) and b{t,x) ~ Hg-7»'(t,2:,i). ^.j^g assumptions placed on 
A and the C"([0,r] x M)-property of u{-, •, 1) imply that Proposition 13.31 can be applied. Hence, 
(f24l) admits a unique solution w(^) G C2+"([0,T] x M x {0, 1}). 

Having established existence and regularity of the solution u of (|2.4p . we turn to the second 
statement - namely, that the process tt*^'^-', defined in the statement, is optimal for (12. 2|) . First, we 
consider the function w : [0, T] x R x M x {0, 1} — ?► (—00, 0), given by 

i;(i,e,x,n) = -e-^(«+"'''(*^^'")). 

The reader will easily check that u is a classical solution of the following PDDE - the formal HJB 
equation for the utility-maximization problem (|2.2p : 

= Vt+ sup [IVxx + TtAwj + ^TT^V^^ + TfVx^ + /iW„ J 

^ (3.21) 

_ v(T,-,-) = _e-^(«+9(-.")), 
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with all regularity inherited from u'-'*'-'. The equation (j3.21|) is the Haniilton-Jacobi-Bellman equa- 
tion for the control problem (|2.2p . where the variables i,,x,n correspond to the wealth process 
/q t^u'^ dSu^\ the Brownian motion B and the jump process iV, respectively. The standard verifi- 
cation procedure (see, e.g., the ideas in the proof of Theorem 8.1, p. 141 in [T5]) can be used to 
show that V is indeed the value function of the utility-maximization problem and that the form of 
the optimal portfolio can be recognized as the optimal value of the parameter tt in the maximiza- 
tion in (|3.21|) . The so-obtained tt is admissible since it is uniformly bounded. The usual difficulty 
one encounters in the course of the verification procedure - namely, the one involved in establishing 
the martingale property of certain local martingales - is easily dealt with here. One simply needs 
to observe that v and ui^'' are uniformly bounded and use the square-integrability of tt e □ 

3.2. Stability of the optimal portfolio. Having established the existence and uniqueness of 
the regular solution to (|2.4p . we turn to the study of the map A i— )■ u'^-' in order to prove the 
statements in the second part of Theorem 12.81 We start with some preliminary growth estimates. 

Proposition 3.5. Let u e C^+^HO^T] X M X {0, 1}) be the unique solution to ()2.4p . Define 

Mo=j\g\o + ^,T\X\l + ^iT, M„ = 7|.g|2+« + (T + l)(i|A|2 + /i). (3.22) 

Then, 

|u(^)|o^Afo, and [m(^)]2+„ ^ M2+"/2e(2+a)Mo^ (3 23) 

Proof. Set Lq = |A|o, La = |A|a and G = I5I2+Q, and let Mq, Ma be as in (|3.22|) . We start with the 
n — 1-slice first; to make the notation more palatable, we write simply u for u('^'(-, •, 1) and A for 
A(-, •, 1) so that u satisfies the linear equation (|3.ip with a = ^^A^ and h = —A. With g = g{-, 1), 
the maximum principle implies that 

\u\o < T\a\o + \g\o < ^TL^ + \g\o ^ Mq. (3.24) 

Corollary IA.5I yields 

[uh+a^G+^[\\ + {l + \\\l+")\u\o 

^G + Ll + il + Ll+'^X^TLl+jG) (3.25) 

^ (1 + l2+")(G + (1 + T)Ll) ^ Mail + i^+") ^ Af2+"/2, 

where the last inequality is a consequence of the fact that the function p (1 -f x''Y^p is nonin- 
creasing in p for a; > and p > 0. Also, due to the first inequality in (jA.Sp and the just established 
(P:^ and ([X^ . we have 

[u]a ^ [u]lf^ \u\F^ < (M'+^)5f^ M^ < m'+^. (3-26) 

We move on to the n — 0-slice, and write u(-,-) for u*^^'(-,-,0) and A for A(-,-,0) so that u 
satisfies the semi-linear Cauchy problem (j3.8|) with a = 2^A^ + l^, h = —A and b = ^exp(— 7-1*). 
Note that u also solves equation (|3.8|) with a = , h = —X and 6 = 0, so that the estimate p.9p 
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implies |e '''"jo = e '''L"-! < e'^l^l" , as well as 

|u|o < max {^T{Ll + 2^i) + \g\o, i log(e^lsl'' + A^r|e-^^|o)^ 

< |5|o + i max (t{^LI + m), log(l + /^T)) . 

This result, together with p.24p and the fact that > log(l + fJ,T), yields the | • |o-bound in 
(|3.23p . By inequality (jA.9p of Lemma IA.6I and estimates p.9p and p.26p , we have 

[e-'^1.<7e^'=Afr^. (3-27) 

Hence, using again the estimate (|3.9p for u, the just established first inequality in p.23p and p.27p . 
we obtain 

So, due to Corollary IA.5[ 

[uh+c ^ [gh+a + ^[X\ + ^[e-^^C^U + (1 + \Ml^"Mo 

^ i(7G + ^Ll + /i) + [e-'^e^^U + (1 + LI+'^Mo 

^ Mo + e^'o+^G^l+l + e''«+^^[uU + (1 + LI+")\u\o. 
Therefore, thanks to Corollary I A. 2 1 and the first inequality in (I3.23p . 

[uh+c ^ e^^+*^« + (e(^G+Af„)(i+f ) ^ ^ ^ ^2+o^ |^|^_ 

Thus, using once more the decrease of the function p i— > {I+x^'Y^p, x, p > and the first inequality 
of (USSI, we get 

The last inequality, along with p.25p , completes the proof. □ 

The multiplicatiye interpolation inequalities of Theorem IA.3I yield the following corollary: 

Corollary 3.6. Let u(^) e C2+"([0,T] x R x {0, 1}) be the unique solution to (HU, and let Mq 
and Ma be as in p.22p . Then, the following estimates hold 

M^']a ^M'+"/'e"*^°, |4-^)|o^Af2+°/2eA/o (3 28) 

We continue with a stability estimate for the nonlinear, n = 1-slice. 

Proposition 3.7. For g e C^+'^{R), e C"([0,T] x K), /i^ e C"([0,T] x M) and 6^=) G 
C"([0,r] X K), /e< uC'), /c 1,2 6e f/ie ttmgue C2+"([0,r] x R)-solution to 

" = + + + a« - 6^6^"'" on [0,T) x M, 

' (3.29) 
^u('\T,-)^g, 

and let 
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Then, 



and 



(3.30) 
(3.31) 



where a(*) = a^^) _ ^ ;j(2) _ /^(i) = 5(2) _ ^(i). 



Proof. Set M*^'') = u^"^^ — u^^^ and subtract the equations p.29p with k ~2 and fc = 1, respectively, 
to obtain that u^'^^ is a C2+"([0,T] x M)-solution to the following semilinear Cauchy problem 



_ 5('5)e7«<" _ 6(i)(eT"<" - e''"*"), on [0,r) x 



(3.32) 



To alleviate the notation, we set 



First, we establish the | • jg-bound in p.30p . Let U{t) = supj.gj{ |u^'*^(t,x)| so that, by the mean 
value theorem, we have 

^W^\t,x) _ ^W^\t,x) < pu{t), for all {t,x) e [0,T] x M. 

Hence, the maximum principle with equation p.32p implies 

u'^^\t,x) < r {Ko + \b^^'>\oPU{s))ds. (3.33) 



J (A'o + |6(i)|oPC/(s))ds. 

We can, therefore, employ Gronwall's inequality (see, e.g., [27j, p. 543) to obtain 

U{t) < (T-Oi^oe'^-*)!"'''!"^ <TXoe^l'''''l"^. 

Noting that Ko < D\h''^'>\o + \a^^^o + P|6(*)|o, we conclude that (jX^ holds. 

Turning to the [•]2+Q-bound (j3.3H) . we use CoroUarv I A. 5 1 and Lemma [A. 61 to conclude that 

^K^ + |6«|oPM'']a + P\u^^^\a\u^'^\o + [5(i)]„P|u(^)|o + D\u^''>\o. 
Therefore, by CoroUarv lA. 21 we have 

[u^%+c. ^K^+ (P|6(^)|o)^+"/' + P\u^'K + Pib^'^U + d) \u(X 
Finally, due to LemmaEl < D\h^^^a + [a'^X + PD\b^^^a, which implies (j33T|) . 



(3.34) 



(3.35) 



□ 



In conjunction with the multiplicative interpolation inequalities of Theorem I A. 3 1 our final result, 
Proposition 13.81 yields the remaining statements of Theorem [ 
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Proposition 3.8. Forge C^+"{[0,T]xRx{0,l}) and k ^ 1,2, letu'-'''' e C2+"([0, T] xRx {0, 1}) 
be the unique solution to (|2.4p . corre sponding to X ^ e C"([0,r] x R x {0, 1}). Then, 

A^o = 7I5I0 + fiT+^Tmax{\X^^^\n, lA^^'lo)', and 
M„ = |5|2+a + (l+T)(^+imax(|A(i)U,|A(2)|„)2). 

Proo/. Let Lo = max(|A(i)|o, |A(2)|o), L„ = max((|A(i)U, |A(2)|„), Go = |g|o and G = |.g|2+o. Then, 

Mo = 7G0 + /iT + iTLg, = G + (1 + T)(l + ). 

Just like in the proof of Proposition 13.51 we deal with the n = 1-slice first; again, to make the 

notation more palatable, we write -u^'^^ for - ,1) and A^'^-' for A'^'^^ (•, •, 1) so that u^''^ satisfies the 

linear equation (jg:^ with a^^) = ^(A('=))^ h'^''^ = -A^'^) and foC^) = 0. Also, set u'^^^ = 

Ao = |A(2) - xW\o and A„ = |A(2) - A(i)|„. Then, u^^^ solves the Cauchy problem with 

/i(^) = _A(2) + A(i), flf'') = ^((A(2))2 - (A(i))2) and fof'') 0. Inequality (13:331) and Corollary[3l] 

yield 

\u^'^\o ^ T(ioAo + |4'^loAo) ^ TAo(Lo + Af^+^/^e^^") ^ TAf^+^/^e^^" Aq. (3.36) 
So, by Corolarv lA.Si estimate p.28p and the decrease of the function p h- > (1 + x^Y^p, we have 
[u^%+c. ^ |z2«UA„+A„L„ + (1 + l2+")|uW|o 

^ (Af^+"/2e(i+")*^" + Af„ + Af i+"/2TAf 2+"/2e*^o^ a„ (3.37) 

^ M4+ae(l+a)Mo^^_ 

Consequently, the multiplicative interpolation inequalities of Theorem IA.3I imply 

^ r2/(2+a)^,f2+2ag(l+a)A/„^^^ (3 

We proceed with the n = 0-slice, and write u^''\-,-) for u^''\-,-,0) and A*-*^-* for A*- •'(•,•, 0) so 
that uC^) satisfies the semi-Hnear Cauchy problem 1^ with af'^) = ^(A^''^)^ + ^, /i^^' = -A^'') 
and ^C^) = ^exp(-7M('=)). Set 

D = Af2+20g(l+2a)Mo ^^^^ p ^ + ^TL^ ^ 

SO that, by Proposition 13.51 CoroUarv 13.61 and inequality p.lSp . D and P dominate (in the ^ 
—sense) the D and P of Proposition [X71 with the current choice of b'^^\u'^^\u (2) and 

Setu(*) =w(2)_u(i). Thenu(^) satisfies the equation with u'l) =u(i),/i(2) ^-A(2)^^(a-) ^ 

-A(2)+A(i),a(*') = 2^((A^'^)'-(A('^)'),&^*^ = ^(exp(-7u(2))-exp(-7?2(i)))and6(i) = ^ exp(-7u(i)). 
Inequality (|3.30|) of Proposition [3771 and the lower bound on u*^^^ in ()3.9[) give us 

^ re^l'''^'l«^(^|(A(2))2 _ (A(i))2|o + D\- A(2) + X^'^\, + P\e-^^''' - e"^''''' |o) 
^Te^'^^ 2 °(Aoio + i?Ao + P|e-^" -e-^" lo). 
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On the other hand, thanks to Lemma lA. 61 the lower bound on u^^^ and m^^-* in p.9j) and inequahty 
p.36p . we have 

Combining the last two inequalities we, obtain 

Im^'^Io ^ Te^^^'"''"^^"'" (Lo + + Pe^^TAf^WZgMo)^^ 

1 2 

^ re'^'''='^''°*^'^''"A.f2+2ag(l+2a)Mo(^ ^_ e^jGo+^TLlj.^^^ (3-39) 

1 2 

where A = log(l + ^T) + A/q + ^Te^TGo+i/zTLj^ + (i + 2a)M^ + 27G0 + \TLI. Since ^ < i, for 
all X > and log(l + ^T) < /iT, we have 

A = log(l + fiT) + Mo + ;iTe''^«+*^«-^^ + (1 + 2a)A/o + A/q + 7G0 - 

^ Afo + eT<=''+*^« + (2 + 2a)A/o + 7G0 ^ e^^^^+i. ^^''^^^ 
It remains to estimate the [-lo+a-seminorm ofu'^^\ By (|3:24| . (|3:25l) and Theorem EH we have 

1 + [u^^^]a ^ 1 + (m2+"/2)"/(2+")Mo'/('+"^ ^ 1 + A.fi+" ^ A/1+". (3.41) 

Applying Corollaries IA.5I and 13.61 and using the decrease of the function p 1— > (1 + a;'')^/'', x, p > 0, 
we get 

+ le-^'''" U|e^^'" - e^i^'" U + (1 + ^i+")|2i^')|o 

(3 42) 

On the other hand, by Lemma I A. 6 1 the lower bound on m^^-' and u^^) \^ (|3.9p . and inequalities 
and dSHD, we have 

^ e-^^«(TA^2+3"/2eA^o _^2.2/(2+a)^^2+2ag(l+a)Mo^^^ (3 43) 

^ M^+3"/2e(2+")AfoA„. 

Moreover, due to Lemma [A . 6 1 and Corollarv l3.61 je"'''"* ' \a ^ 2+0/2^(1+0)1/0 ^nd, with the upper 
bound in p.9p and inequalities p.39p and (|3.40p in mind, 

|e-7n'^) _ ^_^„<i) 1^ ^ ^M„ (TA^2+2age-^o + i^^^^ ^ ^^2+o/2gaMo) ^ [^(5)]^) 
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So, 

\e-'^^'"'U\e-'<^"' - e-^^*"U < Ml+3"e^'"''^'e(2+2")Afo + M2+"/2e(2+a)A/o (3 44) 

Combining inequalities (|3.42p , p.43p and p.44p , we get 

Thanks to Corollary I A. 2 1 we obtain 

[u^%+c. ^ Ml+3"e^"'°*'e(3+2")A^oA„ + ((Af2+"/2e(2+")Mo)i+o/2 ^.^ 1+0/2) 

Finally, due to inequality p.39p . we get 

■7 I Q 2(Mn + l) 

^ Ml+^^e'^ A„. 

□ 

Appendix A. Anisotropic Holder Spaces 

A.l. Definitions and notation. Classical (anisotropic) Holder spaces provide a convenient set- 
ting for stability analysis of a class of utility-maximization problems. Here is a short overview of 
the notation and some basic definitions. 

Let C([0, T] X M) be the set of all continuous functions m : [0, T] x M -> R, and let Cb([0, T] x M) 
be a sub-class of C([0, T] x R) containing only bounded functions. Cb{[0, T] x R) is a Banach space 
under the "sup"-norm 

|w|o = sup \u{t,x)\. 

{t,x)elO,T]xR 

In addition to the "vanilla" norm | • |o, we introduce a family of equivalent, weighted, norms 
{I • |(/3) : /3 > 0}, given by 

= sup e"^^'^-*) \uit,x)\ , for ^ > 0. 

{t,x)e[0,T]xR 

Due to the importance of one-sided bounds, we use the following notation 
\u\ = sup u{t,x), \u\ — inf u{t,x), 

(t,a:)G[0,T]xR (t.x)G[0,T]xR 

as well as their section-wise counterparts 

[^(t, •)] =supit(i,x), [it(i, •)] = inf u{t,x). 

The parabolic distance dp between {ti,xi) and (^2,2^2) in [0,r] x R is defined by 

dp[{tl,Xl), {t2,X2)) = \/\tl - <2| + \X1 ~ X2\. 
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For a function u e C([0,r] x R) and a constant a G (0,1], we define its a-H61der constant 

[u]a e [0,00] by 

u{ti,Xi) - u{t2,X2)\ 

(A.l) 



sup 



(A.3) 



(ti.a;i)#(t2,X2)e[0,T]xR rfp(^(ti,Xi), (^2,2:2)) 

The functional | • |ct, given by 

\u\a \u\o + [u]a, (A.2) 

is a norm and it turns the class C"([0, T] x R) of all functions u E Cf,([0, T] x R) for which [u]a < oo 
into a Banach space. 

For fc G N, the space C''^([0, T] X R) contains all functions u e C([0, T] x R) such that the partial 
derivatives g^„g^„ u exist and are continuous on (0, T) x R, for all m, n G Nq such that 2m + n < k. 
For fc e N and a £ (0, 1] we introduce the space C''+"([0,T] x R) by 

C'=+"([0,r] X R) = {u e C''{[0,T] X R) : -^^u admit extensions in C"([0,T] x 1 

for all m, ri G No such that 2m + n < fc}. 

The norm | • \k+a; given by 

l"U+"= E E l^^wlo for«GC'=+"([0,r]xR), 

2m+n— /c 2m+n<A; 

turns C'=+"([0, T] x R) into a Banach space. In particular, we shall have occasion to use the spaces 
Ci+"([0,T] X R) and C2+"([0,T] x R) with norms 

l^li+a = [u^]a + |u|o + \ux\o, and 

\u\2+a = [Ut]a + [Uxx]a + |ut|o + \Uxx\o + \Ux\q + |w|o- 

Analogous constructions can be performed in the isotropic setting, i.e., in our case, for functions 
It of a single variable. For a G (0, 1], in an act of notation overload, we set 

[u\a^ sup — |u|a = Na + sup . 

xi^x2SM. Fi — a;2| kgr 
Then, the H51der space C'^+"(R), fc G Nq, is a linear space consisting of all functions u : R — > R 
which admit fc continuous derivatives and whose fc derivative satisfies [-^u\a < 00. It 

becomes a Banach space when endowed with the norm | -{k+a defined in analogy with its anisotropic 
counterpart. Further information on Holder spaces can be found in a variety of classical treatments 
of parabolic PDEs; for example, the reader might want to consult for an unbounded-domain 
setting similar to ours, or [22^ and [3D] for a more thorough analysis of linear and quasilinear 
parabolic PDEs on bounded domains. 

In addition to the classical function spaces, we shall have occasion to use functions which depend 
on an additional variable n G {0, 1}. For a G (0, 1], and fc G Nq, the space of all such functions for 
which k + a - Holder continuity is required on both n = 0- and n = I- slices, will be denoted by 
C'^+"([0,r] X R X {0, 1}) (or C'=+"(R x {0, 1}) in the isotropic case). The natural (and Banach) 
norm used in those spaces is the maximum of the Holder norms of the n-slices: 

\u\k+a = max(|u(-, 0)|A;+a, |u(-, l)\k+a), 
[u]k+a = max([M(-, 0)]k+a, [u{-, l)]k+a)- 
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Similarly, for u e C{[0,T] x M x {0, 1}), we define 

M =max(ru(-,-,0)l,K-,-,l)l), and [u\ = min(K-, •, 0)J , L?^(-, •, 1)J ). 

A. 2. Some useful results on Holder spaces. When dealing with various constants in the 
statements and proofs of results below, we use the following convention: 

Convention 2. 

(1) The variables fi, 7 and a G (0, 1) are considered "global" and will not change throughout the 
paper. Any function of the global variables (and global variables only) if called a universal 
constant. 

(2) The notation a <h means that there exists a universal constant C > such that a < Cb. 
Such a constant may change from line to line. 

We start with several well-known interpolation results which we rephrase (and minimally adjust) . 

Theorem A.l (Parabolic interpolation - additive form - [21J, Theorem 8.8.1., p. 124.). 

There exists a universal constant C > 0, such that for any e > 0, and u € C^^"([0,r] x M) we 
have 

Ma <eM2+a + Ce-"/2|y|^^ 

|w,|o <eN2+a + Ce-i/(i+"V|o. 

Corollary A. 2. There exists a universal constant C > such that any function u G C^^"([0, T] x 
M) which satisfies the inequality 

M2+a <D + E\u\q + F[u]a + G|u^|o + H[u^]a, 

for some constants D, E, F.G > 0, also satisfies the inequality 



N24 



<C D 



F 



l+a/2 



2+a 



mo 



(A.4) 



Proof. By Theorem lA.li with the choice of £ = g^, we have [u]a < g^[u]2+Q + C(6i^)"/^ |u|o, so 
that 



Similarly, 



G\u,\o<^[uh+o. + CG'+i+^\u\o and H[u,U < l[uh+c. + CH^+''\u\o. 
The estimate (jA.4|) now follows from 

[u\2+a < D + E\u\o + F[u\a + G\Ux\a + H[Ux]a 



<D+\[u]2+c.+C 



□ 
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Theorem A. 3 (Parabolic interpolation - multiplicative form - plj, Exercise 8.8.2.. p. 125.). 

There exists a universal constant C > such that for any u G C^^"([0,T] x R) wc have 

<CM2+o|^i|r": (A.5) 

Theorem A. 4 (A Holder estimate - based on [51], Exercise 9.1.4, p. 139.). There exists a 
universal constant C > such that for any u G C^^"([0,T] x M) 

[u]2+a < C{[ut + \uxx]a + [u{T , ■)]2+a)- (A. 6) 

Corollary A.5 (A Holder estimate with a transport term). There exists a universal constant 
C > such that for any u G C'2+"([0,T] x M) and h G C"([0,r] x M) we have 

M2+a < C{[ut + iu., + /iuj, + [u{T, Ols+a + (l + |«|o) • (A.7) 

Proof. The Holder estimate of Theorem I A . 41 implies that for f = ut + ^Uxx + hux and g ~ u{T, •) 
we have 

[u]2+a dl [gh+a + [/ - hUx]a ^ [gh+a + [f]a + \h\o[Ux] 

and (|A.7|) follows from Corollary [O] □ 

We finish this section with a preparatory result which states that composition with a C^-function 
X e"^^' is a locally Lipschitz mapping on C"([0, T] x R). Even though it is quite likely that such 
a result is well-known, we were unable to locate a reference, and so, for the sake of completeness, 
a proof is provided. 

Lemma A. 6. For 7 > 0, let : C"([0,T] x R) — > C([0,r] x R) he the composition mapping 

E^u{t,x) = e'^"(*'^) forue C"([0,T] x R). 
Then E^u G C"([0,r] x R) and the following hounds hold for all g C"([0,T] x R), 

|i?7^|o <e^^"l, (A.8) 

[E^u]c.<ie^^^\uU (A.9) 

|£;^m(2)-£;^m(i)|o <7i:'|M^^>-M(^)|o, and (A.IO) 

[E^u^^^ - E^u^'% < 7i^([u(2) - + 7[u(2)]„|u(2) - u(i)|o), (A.ll) 
where D = e-'T"*"! ^T"'''! . 

Proof Since < E'^u < sup^ e'^"^*'^^ (|X8|) holds. For (fO|) . we note that, by the intermediate 
value theorem, 

'e7«(^,y) _e^«(*.^) ^ -fe^^ \u{s,y) - uit,x)\ , (A.12) 
for some convex combination ^ of w(t, a;) and u{s, y), so that e''^ < e'^T"! . 
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In order to get the other two bounds, we pick e C"([0,r] X R), {ti,Xi), {t2,X2) G 

[0,r] X R and set d = dp ((ti, xi), (^2, ^2)) and (5w(*=) = m^''H*2, a;2) - u'-'^'Hti^xi) for fc = 1,2. 
Observe that we can prove (jA.lOp in a similar manner as (jA.9|) . Focusing on (lA.lip , we set 

S = gT"'^' (*2,a;2) _ g7«*^' (*2,a:2) _ ^g7"*^' (*! ) _ gT"'^' (*! ) ^ 

and note the elementary equality 

G{b)-G{a):^{b-a) f G' {h%a,b)) dO, 



where h%a,h) = (1 - e)a + 6b, G e G^{R). Then, 

|,5| = / ^g7''''("<"(ti,a;i),«<'>(t2,a;2)) ^5/ _ ^y(l) / ^^-(h" {u('\ti,xi),u('\t2,X2)) 

Jo Jo 

»1 rl 



^0 -'0 



where 



Clearly, g{r], 9) is a convex combination of values of functions m^^' and u'^^ and the inequality 



5h{e) < -m(i)|o holds. Therefore, 

\6\ < d{[u('^ - «(i)]„7rf" + [u^'^M'^ - u^'Hol^d"), 
which directly implies (jA.lip . 



(A.13) 
□ 
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